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ABSTRACT 



i 

Parameter plane techniques are applied to compensation of third 
order control systems by extending the principle to include the quad- 
ratic form of the parameters. 

Equations for constant bandwidth curves on the parameter plane 
are developed. 

Recursive equations for the extension of these methods to higher 
ordered systems are presented. 

A computer program for the computation and plotting of constant 
zeta, omega and bandwidth loci is written. 
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1. Introduction 



Until 1959 when Mitrovic [i] introduced what is now called the 
Coefficient Plane, analysis and compensation of control systems were 
accomplished primarily by frequency response and root locus tech- 
niques. With the advent of the Coefficient Plane there was then a 
method available which allowed the control engineer to determine the 
root locations of a plant as two coefficients of the characteristic 
equation were varied. 

At about the same time E.R. Ross and T.C. barren [2] in 
conjunction with Dr. G.J. Thaler published a numerical method of 
determining single and multiple section cascade compensation while 
maintaining a constant error coefficient. Shortly thereafter in 
I960 C.D. Poliak [ 3 J and Dr. Thaler published yet a different 
numerical method for stating compensator pole locations for a 
given gain. 

In 1964) Siljak ^4^ extended Mitrovic ' s work to apply to pro- 
blems in feich the system variables or parameters appeared as a 
product in one or more coefficients. This paper will further extend 
the Parameter Plane concept to include the quadratic form of the 
coefficients. Specifically the case of: 

IZ a.s k = 0 (1-1) 

k=0 

where: 

\ • + v* + vy 3 + e k/3 * f i/ 2 + ®k ^ 

and: 

t) k ,C k ,< ^k ,e k ,f k ,g k are rea - 1 - an< ^ cons ' fcan 't 
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o< and are real and variable 

These curves which will be plotted on the cx - /• 
Plane will specify two of the roots of Equation (l 
choice of cx and ^ . 



or Parameter 
■1 ) for any 
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2. Concepts of the Parameter Plane 



The Parameter Plane was derived by a study of the character- 
istic equations of control systems where two parameters were allowed 
to vary. Here a presentation of the basic ideas of Siljak [4] will 
be studied. 

Consider a characteristic equation of the form: 

k 



n 

fn V" 

k=0 



= 0 



( 2 - 1 ) 



where the coefficients a,, a„,..., a are real and s is the 

1 7 d 7 7 n 

Laplace Transform variable: 

s = cr + jo) = - 3 60 + j oi/ 1 - 
By forming powers of s it is shown [4] that: 

s k = 60 k T k (- 3) + iff 1-3^ U k ("^ ) 

where : 

T k (-J) = (-1)\(3) 

V-3) - (-D k+ \(3) 

The and T^ are the first and second kind of Chebyshev 
functions and are generated by the relations M* 

T k+1 (3) -23i k (3) + T k _ 1 < 3 ) 



( 2 - 2 ) 



(2-3) 



(2-4) 



= 0 



-23u k ($) +WS) 



(2-5) 



= 0 



with: 



*o<3> - 
t x (3) = 



u ( 3 ) = 0 



V 3 > 

n- 



9. 



0 



-1 



Taking Equation (2-1), substituting Equation (2-3) in it, then 
equating the real and imaginary parts to zero separately and employ- 
ing Equations (2-4) and (2-5) we will obtain the expressions in 



terms of Chebyshev functions: 
n 



Z7 (-l) L a k UJ k U k _ 1 (^ ) = 0 



lc=0 

n 



ZJ (-i)\«J\(:5) = 0 

k=0 k k 



(2-6) 

(2-7) 



Now, assuming as Siljak did, that the a^. coefficients of 



Equation (2-l) are linear functions of e* and ^3 : 



with: 



®k - V* + \ 

B k ,C k ,( ^k constant and real 

c* and variable and real 

we may then substitute Equation (2-8) into Equations (2-6) and 
(2-7) and obtain: 



(2-8) 



B l* +C 1 



7 s 



B 2 CX + C 2 / / 3 



= -D, 



= — D, 



(2-9) 



where : 



Bl ' S (_ 1 )kb k" 1 Vl 

°i - £ 

h - £ 



n 



k=0 



b 2 = u 

d k=0 k h 



C 2 - £ 

= £ (-D V\ 



( 2 - 10 ) 



10 



Now Equations (2-9) ma y be solved simultaneously for ©t 




giving; 



c?< 



C 1 D 2 ~ C 2 I) 1 
B 1 C 2 “ 3 2 C 1 



A B 2 I) l ~ B 1 I) 2 

r m B 1 C 2- B 2 C 1 



( 2 - 11 ) 



Equations (2-11) now give CX and yS in terms of 5 and CO and 
either may be fixed and the other varied from zero to infinity which 
will give a curve on the cx — ^3 plane. This may be repeated for 
as many values as needed giving an entire family of curves on the 
“-/s plane which completely describe the behavior of a pair 
of complex roots. These curves will usually be constant zeta or 
constant omega loci. 

If s = CT is substituted in Equation ( 2— l) and the same 
procedure is followed, curves of constant sigma or values of s 
along the real axis may be obtained. 

This then is the general procedure to be followed in the devel- 
opment which follows for the quadratic case. 
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3. Cascade Compensation of Third Order Systems 



The initial investigation of the cx — quadratic case is made 
by attempting to compensate a specific third order system with two 
identical sections of cascade compensators maintaining a constant 
velocity coefficient. The uncompensated plant, Figure 3-l> is chosen 
so as to afford not only lead, but lag compensation, and with this 
selection, the system gain is set to place the complex roots at the 
stability limit, which also insures the need for double section comp- 
ensation over a large range of zeta and omega. Establishing the 
generalized compensated plant, Figure 3-2, the characteristic 
equation is determined: 



Simplifying this equation and making a substitution of variables: 



z^s^ + (6z ^ + 2pz^)s^ + (5z^ + 12pz^)s^ + (lOpz^ + 

6p^z^ + 20p^)s^ + (5P^z^ + 60p^z)s + 30p^z^ = 0 (3-l) 




(3-2) 





0 



(3-3) 



and when the terms are arranged in terms of the parameters: 




which shows the general coefficient form: 




12 . 



0 



(3-5) 



G 





K 




T 


S(S+l)($+5) 











UNCOMPENSATED SYSTEM 
FIG. 3-1 



— 



1 1 



Gc 



P (SO- 2} a 




30 

•i 


2 a (S + P) 2 , 




S(S^I)(S^5) J 



COMPENSATED SYSTE M 
FIG^ 3-2 



13 . ! 




or: 



ZJ V3 k = O (3-6) 

k=0 

where all of the c^'s are zero. 

Employing Equations (2-6) and (2-7) 





ZJ (-D^h; (3) . o 

k=0 k k 1 


(2-6) 




ZT (- 1 ) V «o\( 3) - o 

k=0 K K 


(2-7) 


Equation (3-4) may be split into a real and imaginary part in 


terms 


of J and 60 


: 






+ cy3 p 2 f 21 + g 1 = o 


( 3-7 ) 1 


* B 22 


+ 1>2 + + /^^ E 22 + ^2 = 6) 


(3-8) 


where : 






B 21 


= 30^o 2 U 1 (^) 


(3-9) 


B 22 


- 30o5 2 U 2 (J) 


(3-10) 


D 1 


- — 60 coU q ( $ ) 


(3-11) 


B 2 


= -60^U 1 ( J) 


(3-12) 


E 11 


= 2 60 4 U 3 (3) - 126J 3 U 2 ( 5) + 10o) 2 U 1 ( | ) 


(3-13) 


E 12 


= 2 60 4 U 4 ( J ) - 12^o 3 U 3 (^ ) + 10cO 2 U 2 (f ) 


(3-14) 



The second number on the double subscripted variables and 
the number of the single subscripted variables indicates the real(l) 
or imaginary(2) part of the equation. 
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P 21 = + 6 ^\(S) - 5 coU o ( 3 ) +30U_ 1 (|) (3-15) 

P 22 =“^ 3u 3 (^) + -5^(1) + 30U o (J) (3-16) 

G 1 " -*> 5 V $ ) + 6 aj4u 3 ( 3 ) - 5o>\i $ ) ( 3-17) 

G 2 = -oi 5 U 5 (|) + 6oAj 4 ($ ) - 5oj 3 U 3 (J) (3-18) 

Further collecting the terms of Equations (3-7) and. (3-8) we 
obtain: 

c* 2 B 21 + c* P 1 + Q x = 0 (3-19) 

<* 2 B 22 + cXP 2 + Q 2 = 0 (3-20) 

where: 

P-L = (ii-2l)v 

P 2 =/3P 2 (3-22) 

Ql =/3 Ell + y6 2 F 21 + G 1 (3-23) 

Q 2 =/3e 12 + fl\ 2 + g 2 (3-24) 



Tie now have obtained a set of two non-linear algebraic equations in 
terms of $ and cO which must be solved simultaneously for the 
correct values of c*. and ^ . To do this, the method with the best 
chance of sucess appears to be Sylvester's Method in which we form 
a set of four equations by taking the original Equations (3-19) and 
(3-20) and forming two more by a multiplication with o{ giving: 

oc 2 B 21 + <* p i + Qi = 0 (3-19) 

OC 2 B 22 +<*P 2 + 02 = 0 (3-20) 
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0 



( 3 - 25 ) 



c* 3 B 2 i + 



cx 3 B 22 + ck 2 P 2 + o< Q 2 = 0 



Now placing these equations in matrix form; 



’° B 21 P 1 Q l _ 






0 B 22 P 2 Q 2 




2 

c< 


B 21 P 1 Q 1 0 




o< 


B 22 P 2 Q 2 0 




1 



= 0 



we may say that if the ck *s are not zero then; 



0 

0 
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21 



B, 



22 



B 21 P 1 
B 22 P 2 



2 

*1 

Qo 



Q n 

q 2 

0 

0 



(3-26) 



( 3 - 27 ) 



( 3 - 28 ) 



Expanding this determinant we obtain; 

“ B 21 Q 2 + B 21 B 22 Q 1 Q 2 + P 1 P 2 B 21 Q 2 “ P 1 Q 2 B 22 + Q 1 Q 2 B 21 B 22 
- Q ? B 22 " Q 1 P 2 B 21 + Q 1 P 1 P 2 B 22 = 0 



( 3 - 29 ) 



Substituting Equations ( 3 - 21 ) through ( 3 - 24 ) in Equation ( 3 - 29 ) we 
achieve a fourth order equation in ^3 ; 
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4 (-f? 



+ 2F 



22 21 



21 * 22 ^ 21^22 T 22~21 “ ' , 22 " 1"22 



+ D, - F__DfB^_ - 



F 21 B 22 “ F 21 B 2 B 21 + D 1 D 2 F 21 B 22^ + 



/6 3 (- : 



2E 12 F 22 B 21 + 2E 11 F 22 B 21 B 22 + 2F 21 E 12 B 21 B 22 + 
B 1 I) 2 E 12 B 21 ” D 1 E 12 B 22 “ 2E 11 F 21 B 22 “ 



A 



D 2 E 11 B 21 + D 1 I) 2 E 11 B 22^ + 

' E 12 B 21 “ 2F 22 G 2?21 + 2E 11 E 12 B 21 B 22 + 2F 21 G 2 B 21 B 22 + 
2G 1 F 22 B 21 B 22 + D 1 D 2 G 2 B 21 “ B 1 G 2 B 22 “ E 11 B 22 “ 

2F 21 G 1 B 22 - B 2 G 1 B 21 + Wl B 2 2 ) + 



/3(- 2E 12 G 2 B 21 + 2E 11 G 2 B 21 B 22 + 2G 1 E 12 B 21 B 22 " 2E 11 G 1 B 22> + 



(- g 2 B 2 i + 2G i G 2 B 21 B 22 “ G 1 B 22^ 



(3-30) 



from which the coefficients may be determined by a substitution of 
•Eq.ua^iohs (3—9) through (3-18) and the values of the first kind of 
Chebyshev functions, Equation (2-5), in terms of 3 and oh . Since 
the solution of a fourth order equation is at best difficult, it is 
at this point a digital computer becomes a necessity. 

The major problem is not the actual solution of the quartic 
itself, but rather the proper choice of one of the four solutions. 
There are two marked characteristics, however, which help in the 
selection. These are: 

a. Complex answers to the quartic have no physical significance 
and may therefore be discarded as erroneous. 
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b. The definition of cA requires that and ^3 be of 

the same sign so that p and z will be of identical sign. 
Using this information and that available from the Ross-Warren [ 2 ] 
method (Appendix I) as to compensator pole and zero location, it is 
found that the solution to the ^3 quartic (Appendix II) is the 
largest, positive, real value. 

Row entering Euqation (3-21) with this value, and evaluating the 
other coefficients (Appendix II) we find that: 



CX = 



[" Q 1 / B 2l] 



if 



(3-31) 



for in the third order case is always identically zero. 

Thus, with the programing of the appropriate equations, the 
digital computer could give all of the values and plot the constant 
zeta and constant omega loci on the Parameter Plane for any desired 
values. With the <x and ^3 values, Figure 3-3, and the graphs of 
the Parameter Plane Lag, Figure 3-4 > and Lead, Figure 3-5 > compen- 
sation zones, a spot check of the velocity coefficient is made with: 

.2 



K = sG (s) 

v s-^0 c' ' 



= K 




TT z j 
TT p i 



(3-32) 



where K is the compensated root locus gain, 
c 



Case 1. 



= .5 cO = 1.0 



r = -.5 + J 0.86603 



p - /3 « 0.03037 a = p / z = 0.36895 = 0.08231 

c / c c c 



(1.0) (.5* + .86603^) a (4.5^ + ,86603 2 )*( .46963* + .86603^) 
(.41769 2 + .86603 2 ) 



K = 



K = 



*v- 



4.81100 

(4.81100) ( .08231) 2 
(1.0)(5.0)(.03037) 2 



7.06 
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THE ALFA VALUES ARE 



ZETA .000000 .100000 

.50 .23818E-00 .25986E-00 
1.00 .43121E+00 . 45283E*00 


.200000 .300000 
. 27452E*0Q_. 283006*00 
. 45547E*0 0 .441306*00 


.400000 
. 285846*0 0_ 
.411266*00 


1.50 .64650E-00 .65840E*00 
2.00 ,88197Et00 .85706E*00 

2.50 ,11372E*01 .13493E-01 


.62918E*00 .58083E*00 
_,.80 0 0 4E* 0 0_,_72553E*0 0_ 
. 22198E+01 . 473676*01 


.51866E*00~ 

6.33756*00, 

.10206E-02 


3.00 . 14126E*01 . 15121E*01 
3.50 . 17085E*01 .17926E+01 


.173416*01 .209556*01 
_.l 9 4 866*0 1 .21 789 6*0 %_ 


. 258166*01 
. 247296* 0 1 


4.00 . 20252E 4> 01 .21052E+Q1 


. 223736*01 .242116*01 


. 264796*01 


4.50 .23627E-01 .24429E-01 
5.00 .27210E-01 .28037E*01 


.25 6 40E ■ * 01 . 2 7247 E*dl _ 
. 291996*01 , 306796*01 


7291756* 0 1 
. 324106+01 




THE BETA VALUES ARE 






ZETA .000000 .100G00 


.200000 .300000 


. 400000 


.50 . 10 0 89E*0 0 . 865666-0 1 


.738746-01 .623976-01 


.51821E-01 


1.00 .26768E+00 .186896+QO 


. 11578E*00 .50779E-01 


.454986-02 


1.50 . 47270E+0 0 .171566*00 


.17741E-01 .70824E-01 


. 124566*00 


2.00 . 710626*0 0 .51491E-01 


.139976*00 . 227346*00 


. 315036+00 


2.50 ,97780E*00 .477826*01 


. 12213E-02 .27194E*02 


.519356+02 


3.00 .12708E*01 280006*01 


_. 45954E*0l_ . 659l3E*01_ 


.860916+01. 


3.50 . 158626*01 .282146*01 


. 41652E+0 1 . 557606*01 


.697766+01 


4.00 . 192106*01 . 308446*01 


.431246*01 .557576*01 


. 682686*01 


4.50 . 22725E*01_i_34354E*01 


.46 439 JE ♦ 0 1_ . 5 8 7 5 3 E *_0 1 


. 7 09496+0 1 


5.00 . 26384E*0 1 .383336*01 


. 506386*01 .631136*01 


.754826+01 








r io o o 
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.500000 .600000 


.700000 


.800000 


.900000 


2 8 33 4 E +.0 0 275 5 8 EjtA0_. 2 6 2 A3 &♦ 0 0_ 


_.24341 E_±_OQ . 21 753 E±JLO_ 


.36895E+00 .3l500E*00 


.24637E+00 


. 15610E+00 


. 22947E-01 


.44183E+00 . 34833E + 0 0 


. 23376E+0 0 


.87524E-01 


. 12428E+00 


_.5_2_4 2 5 E + QO .39 546EjtO 0_ 


^2 43 57Et 0 0 .58 85 4 E - 0 1 5 4 5 8 8 E ♦ 0 1_ 


.1462QE+Q2 . 12050E+02 


. 79440E+0 1 


. 50060E+Q1 


. 29725E + 0 1 


. 31185E+01 . 35562E+0 l 


. 36999E* 01 


. 33953E-*01 


.25974E+01 


-jt_2_79_8 8 Ejt0.X_._3 0 932 EjtO l_ 


_.3 25 3 2 E 1 01 .3136 6 E* 0 1_ 


_.-2 5 6 6 1 E »Q1 


. 28963E+Q1 . 31247E+01 


. 326Q9E+Q1 


. 31848E + 01 


. 26888E+01 


. 3X255E+Q1 .33167E+01 


.34338E"-01 


.33730E+01 


. 2921 4E + 0 1 



. 34250E+0 1 ,_35934 E j*0 1_._3 6 9 836^0 1 .3 6495E+0 1_._3 2 4 5 3 E+.0 



.500000 .600000 


.200000 


.800000 


.900000 


. 41891E- 0 1 . 32383E-01 


.23076E-01 


. 13717E-01 


.39634E-02 


.30370E-01 .57980E-01 


.89149E-01 


.12779E+00 


. 18641E+00 


^_18 046EjtO 0_.24 0 9 3E*0 0 


31053E + 00. 


9 8 2 E ♦ Q 0_ 


.54396E+0 0_ 


. 40506E+00 . 5 0 0 66E*0 0 


. 60795E + 00 


. 74093E+00 


. 11957E+02 


. 62975E+02 .43894E+02 

| 


. 24815E+Q2 


, 13672E+02 


. 72851E-01 


.1 0318 EjtO 2 ._1 1267E+02 


110 56E-t0 2 *_95651E«-01_ 


.2.0 1 5 2 E + 0 l_ 


.82452E-01 .91940E+01 

1 


. 95796E+Q1 


. 91176E+01 


. 74993E + 01 


.79904E+01 .89518E"-01 


.95378E+01 


. 94838E+01 


. 83371E+01 


.82489 EjtQl_. 92543E + 01. 


• 997 8 IEjLO 1_ 


_^10X9lE^02_ 


.94146E-01_ 


. 87333E+01 . 98020E + 01 


.10649E+02 


.11081E+02 


. 10659E + 02 
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Case 2. S = .2 cO = 3.0 r = -.6 +_ j 2.93943 

/3 = P c = 4.5954 ©< = p c / z c = 1.7341 z Q = 2.6^002 

K = (3.0K.4 2 + 2.93943 2 )~ g (4.4 2 + 2.93943 2 ) If (3.9954 2 + 2.93943 te ) 

c (2.05002 2 + 2.93943 2 ) 

K = 90.21638 

c 

K = (?0,21638)(2.6 ? 002) 2 , 6i00021 

v (1.0) (5.0) (4.5954) 

Case 3. 3 = .8 ^ = 4.5 r = -3.6 Hr j 2.7 

/3 = p c = 10.191 o< = p c / z c = 3.373 z c = 3.02134 

K (4.5)(2.6 2 ± 2.7 2 Al.4 2 ± 2.7 2 )^(6.591 2 +■ 2.7 2 ) 
c (.57866 2 + 2.7 2 ) 

K = 341.3239 

C 

K ,..(M1-32^0..02,1,3.41 2 . . 6>00015 
(i.o)(5.o)(io.19|) 

Uoting that all of the calculated K v ‘s are greater than or equal to 
6.0, the velocity coefficient of the uncompensated system, a check 
of the root loci is made. As shown in Figures 3-6, 3-7 and 3-8, 
the desired root location is also achieved. 
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4. Generalization of the Third Order System 

In Section 3, the development of a specific type one, third 
order system was shown. Here a generalized development of the 
equations necessary to write a computer program for any third 
order astern is shown. 

Starting as before by specifying the compensated plant, (this 
time however in generalized terms), as in Figure 4-1, the character- 
istic equation is written. 



s^ + (2p + 22 r. ) s ^ + (p^ + 2pSr. + 2 TT r. ) s^ + 

i l 2 l 

(p 2 ^^ + 2p2 f jjr i +TTr i )s 2 + (p 2 2, , n , r i + 2p'TTr i )s + 
P 2 Tt ± + K (^2 s 2 + 2-g^ s + p 2 ) (4-1) 



where s 



= the sun of the roots 

= the sum of the product of the roots taken 
n at a time 

TTr. i the product of the roots 
With the substitution of s 




/■? = p cx = p/z (4-2) 

and the collection of like terms in c< and , we obtain: 

c* 2 (Ks 2 ) + c*y^(2Ks) + +2r^s 2 +2^^® + TTi‘ i + K) 

+ /3(2s^ + 22r.s^ + 22TTr.s 2 + 2TTr.s) 

/ 1 £ i i 

+ ( s ^ + 2jt. s ^ +Zj < n i r. s ^ + TTr. s *) (4-3) 

1 g 1 X 
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Employing Equations (2-6) and (2-7) 

B 21 = ) 

B 22 = K 00^(3) 

D 1 = — 2Kci)U q ( $ ) 

D 2 =-2Eo5U 1 (S) 

E 1X = 2 co 4 U 3 ( 3 ) -22^60^(3) 

- 2TTr i ^U 0 (3) 

E 12 = 2«> 4 U 4 (5) ^Zr^^d) 

- 2lTr i coU 1 ( J) 

*2i =-0)^2^^ +Z r i a)2u i(5) 
+ (ITx. + k)u_ 1 (J) 

F 22 = - oo 3 U 3 ( 3 ) + 2 r^ 2 !!^ 3 ) 

+ (TTr. + K)U o (S) 

G-l = -60 5 U 4 ( 3 ) + 2 r i oi 4 U 3 ( 3 ) 

+ TTr i ^ 2 u 1 (3) 

g 2 = -60 5 u 5 (3 ) +2 r i c ° 4w 4 ( ^ ) 

+ 'lTr i oo 2 U 2 ( ?) 

P 1 =^ D 1 

p 2 =/3 d 2 

Ql = /3E n +/ 3 2 f 21 + G, 
q 2 = /(3 e 12 + f 22 + g 2 



and defining: 


(4-4) 




(4-5) 




(4-6) 




(4-7) 


+ 2Z , H 1 r i ^ 2 U 1 (3) 


(4-8) 


+ 2Zpr it 0 2 U 2 (3) 


(4-9) 


-2f^ U o(5) 


(4-10) 


-ZTTr.ooU^S) 


(4-11) 


-S , p n r i <o 3 u 2 ( 3 ) 


(4-12) 


~H r ^ 1 T i U> 3 U i ('S ) . 


(4-13) 




(4-14) 




(4-15) 




(4-16) 




(4-17) 
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we finds 




0 



(4-18) 




(4-19) 



which are two non-linear algebraic equations completely general- 
ized in terms of the uncompensated system poles and^root locus gain, 



(3-30) and (3-3l) are now employed to obtain the solutions for c*. 



Several third order systems were investigated by the application 
of the generalized equations and the Parameter Plane curves, Figures 



the most interesting. Further investigation of three of the curves 
in this family, Figures 4-2, 4-3 and 4-4 shows that there is a rela- 
tionship between K, the root locus gain, o< and jQ . 

These relations are: 



b. Zeta reads directly. 

c. Determine the actual omega at that point by multiplying 
the value read by the cube root of the unconpensa,ted 
system gain. 

d. .Read the value of c< directly from the point chosen. 



S , la and the first kind of Chebyshev Functions. Equations 



and 



, using the same selection rules as before. 



4-2 through 4-9? were plotted. Of these, the K/s^ family appears 



a. Choose a point on the l/s^ 





from the point chosen. 



f. Obtain 




by multiplying this value by 



the cube root of the uncompensated system gain. 



By this method, the values of o{. 




may be determined for 
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all K/s^ systems from the one universal curve. 

Although the majority of the curves are easily interpreted, 
there are some, Figures 4-9 and 4-10, which require careful corre- 
lation of the data points and the curves. As shorn in Figures 4-H 
and 4-12, the correct curves may be picked out and re dr atm if 
necessary for further clarity and ease of use. 
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5. Bandwidth Curves on the Oi - fo Plane 

In many instances, there is also a bandwidth criterion imposed 
on the engineer as well as an optimal operating point for the plant 
under consideration. With this in mind, equations for the plotting 
of constant bandwidth curves on the plane are developed. 

For the purpose of this development a constant bandwidth curve will 
be defined as: 

A constant bandwidth curve for = M is a curve drawn 

upon the parameter plane which specifies the relation between 
the parameters necessary if the transfer function G(s), which 
is a function of the parameters, is to have magnitude M at 
the real frequency 60 ^ . [5] 

Once these curves are obtained they may be superimposed on the 
parameter plane thus indicating what values of the parameters are 
necessary in order to meet the specifications. 

Taking the rational transfer function and defining it: 

m m-1 



G(s) 




p s + p n s 
n ■ L m-1 



+ ...pjS + p c 



n n-1 

V 3 + q n-l s +...1 1 s + q 0 



(5-1) 



where the p ' s and q ' s are of the form: 

n 



p u - Su* 2 + V* * V7 9 + V 3 - y 3 * b 

u — 0,1,2, .. . ,m 

,2 . , ^ .a o2 



% = 



= a^** + b^iC* + c + d ^3 + e ^J3 + f v 



we may write: 



v = 0,1,2, ...,n 

m 



G(s) = 



E 



u=0 

n 



E 4, 



v=0 



(5-2) 



(5-3) 

(5-4) 
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Employing Equation (5-4) we writ© in the parameterized form the 



generalized compensated third order transfer function: 



where : 



and: 







-- 


(5-5) 


P(s) 


2 2 
= «V 


+ 2<a/3ks + /3 2 K 


(5-6) 


Q(s) 


- 


+ Z r fr.s 2 + 2'J*^ i s +ZTT rj + 






/3[2s 4 


+ 22TTr.B3 + 2S7Tr.s 2 + 2 I/DY s ] + 

i 1 £ 1 3 1 J 






[s 5 


+ ETT^S 4 + £TF r iS 3] + H'flVs 2 


(5-7) 



Making the definitions: 



A r = Z] (-l)^ V 6o£a v ; etc. for B r , C r’ D r’ E r’ 



P 



v=0 

even 



v=0 

odd 

m 



s r - i etc. for H ,1 ,J ,K ,1 

u=0 



even 
m 



u=0 

odd 



and substituting in Equation (5-5) > we write: 



G(j^v) = 



(<* <? r + K*) + j(<*/3 I ± ) 



b + / 3 E r + f x ) + i(/3\ + /3e. + F.) 

Setting the magnitude of G(jOi^) = M we find: 



(5-8) 



A i = ZZ (“ 1 ) i '^ V ~ 1 ^ 60 b a v J etc# for B i ,C i ,:D i’ E i’ P i (5-9) 



(5-io) 



0. = £ (-l)^ (u “ l) ^^g u 5 etc. for (5-H) 



(5-12) 
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. I«vl’ . <*> • I - ,i • 'f 



(/3 \ +/3E r + F,) 2 + (/ 3 \ F i ) 

Manipulating Equation (5-13) algebraically we find that: 



(5-13) 

2 



<f> (<*,y3) - H 2 0(<*y3) = 0 



where : 



<f>(<*>/3) = <* V + 2(X 2 K r G r + + <*/<3 2 I? 



(5-14) 



(5-15) 



© (<*y£) - /3 V + 2^ J D r E r + 2^2)^ +y3 
2/3y r + ^ + /3 4 d? + /S^ 2 + ^ + 

2y6 3 E i D. + 2yS 2 D i P i + 2fi\T? ± ( 5 - 16 ) 

Substituting Equations (5-15) and (5-1 6) in Equation (5-14) and 
defining: 

2 2 
P.. = IT + D7 

1 r 1 



E + 2>8"D F + /3 2 E 2 + 



% 

=1 



2D E + 2E.D. 
r r 11 



= 2D F + E 2 + E? + 2D.F. 
r r r i 1 1 



2 2 
- « x i 



V n = 2E P + 2E.P. 
1 r r 11 



W 0 = o< 4 G 2 + 2<X 2 K G + K 2 

2 r r r r 



(5-17) 

(5-18) 

(5-19) 

( 5 - 20 ) 

(5-21) 

( 5 - 22 ) 

(5-23) 



we can write: 
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(5-24) 



M 2 P^ + M 2 Q^3 3 + (M 2 R 1 - R 2 )y(3 2 + + 

(M 2 ¥ x - W 2 ) = 0 

Since the Parameter Plane for compensation purposes has alreacly- 
been determined it is now a matter of taking the computed o< values 
and substituting them along with a constant value of omega and M 
into Equation (5-24) and then solving the quartic. This has as 
its solution the largest, real and positive value of the four roots 
as before. 

Examples of these curves, drawn so as to enable their use as 
overlays on Figures 4-2 through 4-10 are shown in Figures 5—1 
through 5-9 • 
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6. Programing of the Compensation and Bandwidth Curves 



Appendices III and IV are listings of the programs developed 
in the course of investigation of this work. Only PROGRAM PROJECT 
(Appendix III) will he explained, as PROGRAM PLOT (Appendix IV) is 
essentially PROGRAM ROOT LOCUS [6] with a minor subhead and internal 
modification. 

PROGRAM PROJECT is designed to solve the o< quadratic and 
quartic established in Sections 4 and The program is divided 
into two main sections, the first for the computation of the o< 
points and the second for the bandwidth points. 

The first section computes an 80 by 80 matrix of the d and 
points corresponding to set values of 5 and <■& . The computa- 
tional part is followed by two distinct graphing sections, one for 
lag and the other for lead compensation. 

The lag graphing section is set up so that during the plotting 
of the curves each value of e< is tested to determine if its value 

H 

is 10 £ ^ 1.0001 . If no points are found within this range 

then a print out is made: 

MO LAG COMPENSATION POSSIBLE 



For the lead section graphs, is again tested by the criterion 

1.0001^ ^ (X-graph scale) (X graph width) . Again if there are no 

values of o( within this region the statement: 

M0 LEAD COMPENSATION POSSIBLE 

is printed. In this case however a study of the printed values of 
©( must be made to insure that the points are indeed non-existant 
or rather just lie outside the range of the graph. 

The second main section of the program computes the value of 
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/3 for a given value of «, as determined by the X graph scale. 
Here the plotting routine is set up so as to not plot zero points 
and to stop the curve when either the o< or value exceeds the 
range of the graph. 

There are several go-no go options available within the 
program which are amply explained in the comment section preceeding 
the program listing in Appendix III. 
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7. Extensions to Higher Order Systems 

Although the work presented to this point has been limited to 

third order systems and the program written for this specific case, 

investigation shoirs that generalized equations may be written which 

will allow the extension of the program to higher ordered systems. 

't/ll 

It can be shown that for a given n order system with no zeros to 
be compensated with two identical sections of cascade compensation, 
that the characteristic equation of the system may bo generally 
written as: 

n+2 0 n+1 2 n / 2 2 — 2\ Tr 

s + 2ps +ps + (z s + 2pzs + p )K + 



J=n J=n 

k=l _ . „ k=0 

2p T, (ETPr.^ + p 2 ZZ 

k=n 3 k=n-l 

0=1 0=1 



(ETTrjJ 
a 



k 



j=n 

k=2 



zr (Zt/vk 



= 0 



(7-1) 



k=n+l 

0=1 



where for n=4 the equation would be written: 

s^ + 2ps^ + p 2 s 4 + (z 2 s 2 + 2pzs + p 2 )K + 

SpCETfr.s 4 + ^TTr-s 3 +2^ r.s 2 +ZV' ? ± s ) + 

p 2 (Z'ITr i s 3 +£ , n l r i s 2 + 2 , n 1 r i s + 2TTV) + 

(STFr.s 5 +£, r n 1 r i s 4 +£ r jr r i s3 + Zfr ± s 2 ) = 0 (7-2) 

It may be further shown that the parameters defined by Equations 
(4-4) through (4-13) may be written: 
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(7-3) 



b 21 = Ko> 2 U 1 (3 ) 



B. 



22 



Kco^UgCS) 



\ = -2ICo5U o ( $ ) 
D 2 * ^KajU^ $ ) 



(7-4) 

(7-5) 

(7-6) 



0=n 

k=l 



E. 



11 



2(-l) n+1 .o n+1 U n (5) + 2 IT [(-l) k a;V 1 (S)(ZTPr i )1 

k=n u - J 

0=1 (7-7) 



0=n 
k=i r 



E. 



12 



- 2(-D n+1 « J n+ \ + i< 3 > + 2 c[(- 1 ) k «\(S)(2'P’ r i)] 



k=n 

0=1 



(7-8) 



P 21 * < 



0=n 

k=0 



-l) n 40 n U n _ 1 ( S ) + E [(-!) k 60 k U k _i( 5 ) (2 7Tr i )J + KU_ 1 ( 5 ) 

^ (7-9) 



0=1 



0=n 

k=0 



F, 



22 



(-l^cO^CS) + £ r (- l ) k O )\(^)( r ,? r i )l +KU o (S) 

k=n-l u J J 



0=1 



(7-10) 



0=n 

k=2 



- (-D n+2 <o” + \ +1 (S) ♦ ,2, 



k=n+l 

0=1 



(7-H) 



G 0 = ( 



0=n 

k=2 



-D n+2 a3 n+ \ +2 (S) + c [(-i) fc o 4 \(S)(ET» i j] 



k=n+l 

0=1 



(7-12) 
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These then are the recursive equations required for the complete 
generalization to a n ^ 1 order system. By employing the above 
equations and replacing in PROGRAM PROJECT (Appendix III) cards 100 
through 150 and 300 through 540 with the appropriate programing, the 
program may be used for any given n^ k order system. 

In like manner by generally defining: 

. 2 ,, 2 



and: 



P(s) = o< 2 Ks 2 + oy3 



Ks 



y3 2 K 



„/ \ n+2 

Q(s) = s + 



0 =n 

k=l 



2/S s n+1 + /S 2 s n + 2 fZ (STTr.)s k + 
' ' k=n i 1 



0=1 



(7-13) 



0 =n 

k =0 



0=n 

k =2 



s 2 in (E?r.)s k + ^ (2 f lT r i) s ' 

k=n-l 3 k=n+l 3 



0=1 



0=1 



(7-14) 



and using Equations (5-8) through (5-11) we may replace in the 
program cards 2860 and 2880 through 2920 , thus adapting this part 
of the program to a general n order application. 
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8. Conclusions 

Throughout this development of the Parameter Plane quadratic 
extension, the c^'s in the generalized coefficient form: 

(b k o( 2 + c k c< + djoyg + e^3 + f^/3 2 + g k ) =0 (8-1) 



have been identically zero. This at first appearance might seem to 
detract from the generalization. The inclusion of this parameter 
does not however introduce any great difficulty in the solution. 

The change in the development would be to the value of P^ and P^ 
which would become: 



and the final 

CK 



h ■ c i + M 

P 2 - °2 + P\ 

solution for c* 




which would change to: 



P I “ 2B 21 Q 1 



X 

3 



4B 



21 



( 8 - 2 ) 

(8-3) 



(8-4) 



For this case, new selection rules for acceptable values of ©< 
would be used, and would be much like those presented for jQ . 

Though the extension of the Parameter Plane to include the 
« - quadratic case makes this tool even more useful, further 
work is still to be done in this field. Not only must the equations 
for the solutions of the Parameter Plane curves for such cases as: 



- V 2 ^ 2 + <y*/^ + + v< 2 + fy3 2 * 

e k <*/3 + + i^/3 * \ [5] (8-5) 



60 



and higher ordered combinations of the parameters be developed, but 
more efficient programing techniques must be developed. In the use 
of PROGRAM PROJECT, for instance, as the location of the system poles 
on the O” axis of the S-plane move to the left, the computational 
time becomes excessive due to present programing technique and computer 



Another major problem in further extensions of these techniques, 
and indeed even other applications of the curves from the preceeding 
development, will be interpretation. In this case, the initial sub- 
stitution of variables immediately allowed interpretation of the 
curves sight unseen. Here then, will be most likely the one single 



used as representations of other variables in control systems, each 
application will have its own unique interpretation. 



speed. 



drawback to further extension, for as the parameters <=* 




are 
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APPEHDIX I 



THE ROSS-WAHEEN METHOD OP COMPENSATION 
The Koss-Warren C2] Method is an algebraic way of determining 
the location of compensator poles and zeros for single or multiple 
sections of cascade lead or lag compensation. The relationships 
of the method are based on relationships between the root locus gain 
and the error coefficients, where the uncompensated root locus gain 

is ’ , , M n t> + piI 

| K «| • m f 1 " 1 ) 

TT |r + z \ 

1=1 

and the compensated root locus gain at the point r is: 



K 0 (r) 




( 1 - 2 ) 



where k is the number of identical sections introduced by the 
compensator. 

The other part of the equality for the derivation is determined 

by the fact that the uncompensated error coefficient is defined by: 

m 



K x (r) 



K(r) 



r 

i=l 



n 



(1-3) 



TP Pi 

i=l 1 



and the compensated error coefficient is: 



K (r) = K (r) • 
ex'' y x' ' z . 



(1-4) 



d ^c 

These relationships are expressed geometrically on the S-plane in 
Figure 1-1 and are used to derive the relationships: 
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s- PLANE RELATIONSHIPS TO 
DERIVE THE TH£ ; ROSS- WARREN EQUATIONS 
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hit sin 

z b “ k(r) sin 0 



(1-5) 



where: 



and: 



Here: 



_ h sin A 
p c " sin (© - <J>) 



K 



= cot”" (cot $ " K^y CSC 

© - 'Tf + <j) - ^ ^ 



4> 



(1-6) 

(1-7) 

(1-8) 



h - the distance from the origin to the root point 

and is specified by the choice of zeta and omega 

K - root locus gain of the uncompensated system 

K(r) - root locus gain of the uncompensated system at 
the point r 

<|) - angle required to make r a point on the root locus 

^ - an angle calculable from K, K(r) and <|> 

3 - angle between the negative real axis and h 
© - an angle calculable from <J) , (f and 



These equations are generalized for n sections so that: 




(1-9) 



( 1 - 10 ) 



where : 


- 


cot ^ (cot <|> n - 


yicTrT 080 


(1-11) 




e„ - 


« * 4n- s 


- >n 


(1-12) 


and: 




II 

P 




(1-13) 



65 



APPEIIDIX II 



DETERMINATION OP BETA QUARTIC SOLUTION 
Here, throe examples of the approximately fifty calculations 
that were made in the investigation of the proper solutions of o< 
and , will he presented. Those calculations are started by em- 
ploying Equations (l-9) through (l-13) with n equal to two* 

Case 1* 3 = 90° K = 30.0 h = 4*0 

4> = tan -1 4«P . tan -1 ^ + tan -1 ^ - 180° = 24° 38' 

T 1#U 2#\J 

<j> 2 = 12°19» 

.1 A 

K(r) = (4.0) ( 17*0) 2 (41*0) 2 = 105*60304 

YK/K(r)' = *532994 



>\ ? = cot _1 (cot 12°19’ - (. 532994) esc 12 u 19’) 



25°40‘ 



0 = 180° + 12°19 ’ - 90° - 25°40‘ 



76°39 i 



_ 4.0 sin 25 40* 
sin 64°20« 



= 1.9225 



(4.0) (>532994) = .94906 



sin 76 39 » 



this implies that: 

/3 = 1.9225 c< = 2.02568 

Case 2. ^ = 60° K« 30.0 h = 5.0 

<f> . tan -l(^|22) . tan (-4^02) . tan (jUi°R) . l8o ° 

= 109°16» 

<j> 2 = 54° 33' 

K(r) = (5.0) (21.000199)^(25.00019)“ = 114.8764 
^K/K(rj^ = .511029 
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7\ 2 = cot -1 (cot 54°33* - (.511029)csc 54°33') = 85°09.5' 
© 2 = 186° + 54°33' - 60° - 85°09.5' = 89°32.5* 



(5.0) sin 85°09»5 t 
sin 34°50.5' 



8.72057 



z = (5.0) (.511029) = 2.54616 

sin 89 23.5’ 



this implies that: 

/3 = 8.72057 = 3.42498 

Case 3. 3 = 90° K = 30.0 h = 5.0 

(|> = tan ^(q^q - ) + ‘t an ^(^q) + ‘tan 1 (^§) — 180° = 33°48* 

4? 2 = 16°54' 

1 1 _ 

K(r) = (5.0)(26.0) s (50.0) s = 180.277586 

= .407934 

)\ 2 = cot -1 (cot 16°54' - ( .407934) esc 16°54') = 27°59' 

0 2 = 180° + 16°54» - 90° - 27°59 * = 78°55’ 



(5.0) sin 27°59* 



P c " o 

sin 62 01» 



= 2.60005 

z_ = (5.0) (.407934) - s - n 2 ? — = .95444 



sin 78 55* 

this implies that: 

yS = 2.60005 = 2.72416 

Now taking the values of just found we compare them to the 

roots of the y3 quartic given in Figures II-lb, II-2b and II- 3b 
respectively. We see that the proper answer for y3 , in all cases, 
is indeed the largest positive and real root of the quartic. From 
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Figures II-4 and II-5 we now obtain the values necessary to solve 



for the value of . 
Figure II-6, and again 
selected values of ^ 



These steps for the throe cases are shown in 
we note that only when applying the above 
do we obtain the correct value of , 
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VALUES OF i HE BETA COEFFICIENTS FOR ZE'fA = .00000 
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values or THE beta coefficients for ZETA = .50000 
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C***CARD FIFTEEN - BWX.BWY - THE SAME AS CARD EIGHT EXCEPT FOR THE 
C BANDWIDTH CURVES. 

C BWY IS ALSO USED TO DETERMINE WHICH VALUES OF ALFA WILL BE USED IN 
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IF ( N-20 )204»205»205 003280 

205 MODE = 3 , 003290 

204 CONTINUE 003300 

DO 212 I =1 » 20 003310 

IFCBFIN (N* I > - .000001 )212»212»209 003320 
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456 IF(E(NP3-1 ) )458»457,458 00004660 
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B ( 6 ) = BETA*BETA 000300 
CALL PLTRTS (A,B) 000310 
GO TO 3 

END 000320 
SUBROUTINE PLTRTS < A * B ) 000330 
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